D-Wave Checkerboard Order in Cuprates 
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We show that a d-wave ordering in particle-hole channel, dubbed as d-wave checkerboard order, 
possesses important physics that can sufficiently explain the STM results in cuprates. A weak d- 
wave checkerboard order can effectively suppress the coherence peak in the single-particle spectrum 
while leaving the spectrum along nodal direction almost unaffected. Simultaneously, it generates a 
Fermi arc with little dispersion around nodal points at finite temperature that is consistent with 
the results of ARPES experiments in the pseudogap phase. We also show that there is a general 
complementary connection between the d-wave checkerboard order and the pair density wave order. 
Suppressing superconductivity locally or globally through phase fluctuation should induce both 
orders in underdoped cuprates and explain the nodal-antinodal dichotomy observed in ARPES and 
STM experiments. 

PACS numbers: 74.25. Jb, 74.25. Dw, 74.75.-h 



Recently, Scanning Tunneling Microscopy(STM) has 
revealed surprising yet important electronic structures 
in the high temperature superconductors. The fourier 
transform scanning tunneling spectroscopies (FT-STS) 
from STM have captured two different general features 
in both momentum and energy spaces^^^^^^^^^^^^^^^^^. One 
feature is dispersive peaks in FT-STS^i^, interpreted as 
interference patterns caused by elastic scattering of quasi- 
particles from impurities^^. The other is non-dispersive 
peaks, a checkerboard modulation observed in various 
different materials and circumstances. The checkerboard 
structure was first discovered locally in BSCCO near a 
vortex core.^'^^ Then, it was found to be a characteristics 
of the large gap regions where the STM spectrum resem- 
bles that in the pseudogap phase. ^^^^^ Later, the STM 
studies of Ca2-a?Naa^Cu02Cl2 revealed the presence of a 
global long range commensurate checkerboard order in- 
dependent of doping^ Finally, in the pseudogap phase, a 
similar checkerboard pattern was also observed.^ 

The origin of the checkerboard has become central to 
understanding the nature of electronic states in cuprates. 
Various different mechanisms have been considered to ex- 
plain the observed non-dispersive checkerboard modula- 
tion, including pair density modulationi^ii^ii^ii^ii^, cur- 
rent density modulatio n^ ''i^^ , spin modulation^^ , stripe 
charge modulation^^'^, and impurity scattering^^ and 
so on. Among the proposed mechanisms, the pair den- 
sity wave (PDW) has been shown to capture important 
characteristics of the checkerboard density modulation. 
The mechanism of PDW derives from high pairing en- 
ergy scale in cuprates. It suggests that unlike the super- 
conductivity of normal BCS type superconductors that 
can be destroyed by breaking Cooper pairs, the super- 
conductivity in cuprates can be more easily weakened or 
destroyed by phase fluctuations than by pair breaking. 
Based on this argument, pair density localization^^ was 
first proposed to explain the local checkerboard modu- 
lation in the presence of impurity or vortex. Later, a 
global pair density wave (PDW) was proposed to explain 



the checkerboard physics in the pseudogap state^^^^. It 
has also been shown that the symmetry of the tunnel- 
ing intensity can distinguish the pair density modulation 
from the conventional density modulation^^. While the 
pair density modulation provides a good understanding 
of the experimental results, it does not establish a direct 
link between superconducting and pseudogap states, as is 
suggested by the presence of the non-dispersive checker- 
board density modulation in both states. Furthermore, 
the theory has not explained three important features 
present in the STM spectrum: 1. the density of state 
at low energy in superconducting state does not change 
whether or not a checkerboard modulation takes place lo- 
cally; 2. the overall intensity of the modulation is rather 
small; 3. the small modulation has a large effect on the 
STM spectrum around the superconducting gap. 

In this Letter, we show that an explanation based on 
d-wave checkerboard density (DWCB) order in particle- 
hole channel can overcome the above limitations of the 
PDW theory. The DWCB can be viewed as a natural 
extension of the d-density wave(DDW) order proposed 
to explain pseudogap physics^^^^^, and is only differ- 
ent from the latter in terms of order wave vectors. We 
show that the DWCB order must exist when the PDW 
order is present in the global d-wave superconducting 
state. Moreover, we demonstrate that the DWCB, with 
Q = {(7r/2,0),(0,7r/2)} and /(k) = cos(A:^) - cos(A:^), 
can explain the STM experimental results. We show 
that the DWCB order has little effect on the density of 
state at low energy in the superconducting phase, but 
has a strong effect on the STM spectrum around the su- 
perconducting gap at high energy. This result natually 
explains the puzzling dichotomy between the nodal and 
antinodal regions observed in STM'^ and angle resolved 
photoemission spectroscopy (APRES)^. The DWCB or- 
der also preserves in FS-STM spectrum at the wavevector 
Q the same symmetry as that observed in the experi- 
ments. Moreover, the DWCB preserves the nodes in the 
single particle spectrum, and generates a Fermi arc with 




FIG. 1: The configuration of bond density of DWCB order, 
or Ws (r) . It is manifestly shown that the pattern has 4a x 4a 
periodicity and (1^2 _y2 symmetry. 



little dispersion around the nodal points at high tempera- 
ture, which are consistent with the results from ARPES. 
The Fermi arc has been a signature of the pseudogap 
region, and has been proposed to explain the checker- 
board pattern observed in the pseudogap state^^. Thus, 
the DWCB provides a physical origin of the Fermi arc. 
To our knowledge, this letter presents the first concrete 
model describing this physics. 

The connection between particle-particle (P-P) and 
particle-hole (P-H) channel orders: To illustrate the com- 
plementary connection, we make use of the DDW or- 
der, since the only difference between the DWCB and 
DDW orders is the order wavevectors. Let's consider a 
state with a DDW order, (Z)^ cJ^Ck+Q'a) = iAddw/(k), 
and a d-wave superconducting order (DSC), (ck^C-ki) = 
Adsc/(k), where Q' = (7r,7r). Some simple calculation 
will show that in the above mixed state of DDW and 
DSC, there naturally exists a PDW order with a wavevec- 
tor at Q', given by 



(ckTC-k+Qu) ^ ^AddwAdsc/^(k). 



(1) 



This indicates that the mixed state of DDW and DSC has 
a complementary description as a mixed state of PDW 
and DSC. It is important to note that the symmetry of 
the PDW order in this case is an extended 5-wave. It is 
also easy to see that the symmetries of the order in the 
P-P channel and that in the P-H channel must be cor- 
related with each other: if one is the extended 5- wave, 
the other is the (i-wave and vice versa. The above result 
holds for the DWCB order by replacing Q' by Q. In a 
global DSC state, the DWCB order must exist when a 
checkerboard PDW order is present. The complementary 
connection suggests that if the phase fluctuation leads to 
a Cooper pair modulation pattern, orders in both the 
P-P and the P-H channel have to be simultaneously con- 
sidered in microscopic models. 

BDG equation and co- existence of PDW and DWCB: 
To verify the existence of the DWCB order and the above 
complementary connection, we perform a self-consistent 



BDG calculation. We start from the following general 
Hamiltonian on a two dimensional square lattice with 
the nearest-neighbor attractive density interaction, 

H = -i^Y^ [Ujclcj + h.cj -^"^VijUifij - ^Y^m. (2) 
The density interaction Vij includes two parts. 



V^J 



(3) 



where V^ favors a d-wave superconducting state and 
^yvi.Yj describes a modulating density interaction which 
creates a small checkerboard modulation on the top of 
the uniform superconducting state. 



8Vr 



l///_r(cos Q • r + COS Q • r'). 



(4) 



where Vi = V and VI = —V . The difference on the 
signs of V^ and V^^ provides us the d-wave symmetry in 
the P-H channel order. We also note that the sign differ- 
ence does not break rotational symmetry with respect to 
a proper rotation center in the lattice. 

Starting with Eq. ([2]), we can derive the BdG equa- 
tions by introducing mean- field decoupling of the nearest- 
neighbor interaction terms and obtain a self consistent 
solution. The Bogoliubov-de Gennes equation is given 
by 
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(5) 



where Hq and A are transfer matrices such that 

Fo^n(r) = - ^^{t + Ws{v)) + M(r)] Vn(r + S) (6) 
A^„(r) = ^^A,(r)^„(r + ,5) (7) 

where ijjni}^) can be either Un{v) or Vn{Y)^ and 5 denote 
nearest-neighbor vectors. The order parameters are self- 
consistent ly determined by the self consistent equations: 
the (i-wave pairing amplitude on a bond (r, r + ^) is given 

by 



AW(r) = T/0(Cr,c 



r+.5T "•" '-r+Sl 



CrT)/2> 



(8) 



the pair density wave order in the P-P channel is 

,(2) 



A^^^(r) = 5V,,,+s{c,^c,+,^ + Cr+,iCrT)/2, (9) 



and the density order in the P-H channel is 

W5{y) = -SVr^r+5{cl^Cr+5a)' 



(10) 



We have numerically solved the BdG equation with 
various different parameter settings including different 
band structures in the square lattice with different N xN 
sizes. We find that the co-existence of the PDW and 
DWCD and the symmetry correspondence between them 
are the robust results in our calculation for this system. 



For example, with a parameter setting, V^ = 2.5t and 
V^ = t, the results are given by, for r' — r = f or ^, 

A^'LW = ±Ao, A^^IrW = Ai(cosQ • r + cosQ • r^, 
andl^r'-rW = ±iyo(cosQ-r+cosQ-r') with Aq = 0.3t, 
Ai = 0.25t, and Wo = O.lbt , which are corresponding 
to a state with 



(a) 



Adsc(r) 

Apdw(r) 

Wdwcb(r) 



Ao 

Ai cos Q • r 

Wo cos Q • r 



(11) 

(12) 
(13) 



We note that PDW and DWCB has the same spatial 
modulation as cos Q-r, while the symmetries of the PDW 
and DWCB orders are the extended 5-wave and d-wave 
respectively. As shown in Fig. [U DWCB has 4a x 4a 
periodicity and (i-wave symmetry. Similar order param- 
eters have been mentioned in Ref— . The solutions are 
independent of the initial guesses for the local variables 
and converge quickly as A^ increases. 

After demonstrating the co-existence of the DWCB 
and PDW, now we are interested in the physical ef- 
fects of the DWCB order. To obtain a clear picture of 
the DWCB, we illustrated a static pattern of the bond 
strength of the DWCB order in Fig. [T} 




W^{r) = Re{Wo) 
Wy{r) = -Re{Wo 
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FIG. 2: (a) The average DOS in the particle-hole symmet- 
ric case. The inset shows the average DOS with the finite 
chemical potential included in the band dispersion provided 
by Norman et al^^ . Aqp represents Wo- (b) The Fourier 
components of LDOS at Q = {(7r/2,0), (0,7r/2)}. 



as such 



It is clear that the DWCB order has 4a x 4a periodicity 
and dx2-y2 symmetry. Similar order parameters have 
been mentioned in Ref^. 

Now we numerically calculate the average density 
of states(DOS), p(cj), and the Fourier components, 
PQ(a;), at the wavevectors of the DWCB order: Q = 
{(7r/2,0), (0,7r/2)}, and directly compare them with ex- 
perimental results. We calculate the above quantities in 
two situations with different band dispersions. The re- 
sults are rather general and are insensitive to the bare 
band structures. First, we performed calculations in the 
particle- hole symmetric case. For a simple band disper- 
sion, we choose t = — 125meV and /i = 0. Ao = 40meV, 
which is relevant for under doped BSCCO. The imagi- 
nary part of the self energy 77 = 5meV is used for the 
entire numerical calculation. Fig. ^a) shows the aver- 
age DOS(Q = (0,0)) normalized by the non-interacting 
Fermi liquids. In the absence of DWCB, there are sharp 
coherence peaks at the energy of superconducting gap. 
As DWCB order develops, the coherence peaks are sup- 
pressed and pushed away while the spectrum at low en- 
ergy remains unchanged. Fig. El^b) shows the Fourier 
components of the local density of states(LDOS) at the 
wavevectors, Q. As expected, Pq{uj) is even with respect 
to cj, namely, Pq{uj) = pq{—uj). Second, we repeat our 
calculations with the band dispersion provided by Nor- 
man et al?^ and the result is displayed at the inset in 



ti/2(cos kx + cos ky) -\-t2 cos kx cos ky 
ts/2{cos2kx + cos2ky) 
^4/2 (cos 2kx cos ky + cos kx cos 2ky) 
^5 cos 2kx cos 2ky — /i. 



(16) 



where h = -0.5951eV, ^2 = 0.1636eV, ts = -0.0519eV, 
U = -OlllTeV, and t^ = O.OSlOeV. The chemical po- 
tential /i is now set to -0.1660eV. Compared with the 
particle-hole symmetric case, the effect of DWCB on the 
LDOS is insensitive to the energy band structure. Qual- 
itatively the numerical results are strikingly consistent 
with experimental results^^^, and the large gap region 
can be represented by the presence of DWCB which is 
weak at 8-12meV. 

Analytically, the general features in STM measure- 
ments can be captured by the DWCB. First, due to the 
anisotropy inherited from the d-wave factor of pairing, 
a weak DWCB order has a much stronger effect on the 
antinodal region than on the nodal region, thus naturally 
explaining the puzzling dichotomy between the nodal and 
antinodal excitations in high-temperature superconduc- 
tors: The local phase fluctuations of Cooper pairs lead 
to a local modulation of (i-wave ordering in the particle- 
hole channel, which strongly affects the antinodal single 
particle excitations. Second, like PDW, DWCB is bond- 
centered and consequently, Pq{u;) is an even function of 
Fig. [2l^a). The band energy dispersion is now modified co, toc^^. The symmetry has been shown to distinguish 



the PDW from the typical particle-hole CDW. The ex- 
isting experimental results are consistent with the even 
case. 

The above result demonstrates the consistency be- 
tween the DWCB order and the STM experimental re- 
sults in superconducting state. Now we show that the 
DWCB order also captures important physics in the pseu- 
dogap phase. One important feature of the pseudogap 
phase is the non-dispersive Fermi arc developed from the 
nodal point along Fermi surface observed in ARPES^. 
The Fermi arc has been used to explain the STM re- 
sult in the pseudogap phase^^. If the pseudogap phase is 
strongly connected to phase fluctuations of d wave super- 
conductivity, the single particle spectrum should reflect 
the DWCB order. Therefore, a robust Fermi arc feature 
should exist in the mixed DWCB and DSC phases in high 
temperature. We found that this is indeed the case. 

In Fig. [31 we have calculated A(k, co) in the pseudogap 
state within the model of Franz and Millis^^. Fig. [3](a) 
shows the numerical solutions of the spectral weight, 
74(k, cj), as a function of a; at high temperature T = 120K 
where Aq = 40meV and the DWCB order is equal to 
8meV. As expected in the pseudogap state, the scatter- 
ing vector connecting the tips of Fermi arcs unchanges 
as energy uj increases, and is nearly equal to wave vec- 
tors of the ^-checkerboard order parameter, |Q| = 7r/2. 
Fig. [31(c) shows the spectral weight at cj = for the cuts 
perpendicular to the Fermi arc which matches the exper- 
imental results^^. For the temperature dependence, the 
length of the Fermi arc is linearly increasing as the tem- 
perature rises above Tc^^. As seen in Fig. [31(c), at very 
low temperature the Fermi surface is gapped except for 
nodal point, (7r/2,7r/2). As the temperature rises, the 
gapless region elongates along the Fermi surface, with 
slight broadening in the direction perpendicular to the 
Fermi surface. 

In summary, the DWCB order offers a unified explana- 
tion for the STM experiments in both the superconduct- 
ing and pseudogap phases. A number of important issues 
need to be addressed. First, while we show that the com- 
plementary connection exists between the d-wave order 
in the particle-hole channel and the pair density wave or- 
der in the particle-particle channel, the order wavevectors 
should be determined by microscopic models. In gen- 
eral, the pair fluctuation in different mircoscopic models 
can lead to different order wavevectors in the particle- 
hole channel. For example, the pair fluctuation can be 
anisotropic in the space that breaks the rotational lattice 
symmetry, and will result in a stripe-like one dimensional 
ordering. Second, it is interesting that in a recent renor- 
malization group study of the electron phonon interac- 
tion in cuprates, the authors has shown that the DWCB 
order rises from the coupling to the half breathing or the 



Big phonons^^. Combining their results with ours sug- 
gests that the superconducting phase fluctuations may 
be strongly coupled with phonons. This hypothesis is 
of great importance and in need of future investigation. 
Finally, although a full calculation based on the BDG 
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FIG. 3: (a) With DSC and DWCB coexisted, the spectral 
weights, A(k,oo), are plotted as a function of oo (-30meV ^ 
OmeV) in the first Brillouin zone. Here we used t — — 125meV, 
t' ^ II ^ OmeV, and T = 120i^. Each BZ is segmented 
by 100x100. White dotted lines indicate the non-dispersive 
Fermi arcs, (b) Temperature dependence of A(k,a;) at the 
Fermi level, uj — OmeV. (c) EDC curves from the nodal to 
antinodal point along the Fermi surface show the gapless re- 
gion, or Fermi arc. 



equation with the DWCB order is yet to be completed, 
preliminary results show that the qualitative conclusion 
drawn here should remain valid. 

In conclusion, the order in particle hole channel, i.e., 
the DWCB order, can explain the STM spectrum and the 
nodal-antinodal dichotomy observed in both STM and 
ARPES experiments^^^^. The presence of the DWCB or- 
der also preserves the gapless dispersion at nodal points 
and simultaneously generates a Fermi arc with little dis- 
persion around nodal points at finite temperature. The 
results are consistent with ARPES experiments and pro- 
vide an explicit physical explanation of the fermi arc in 
the pseudogap phase. 
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